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ABSTRACT 

We present moderate and high resolution infrared spectroscopy of the peculiar eruptive 
variable V838 Mon, which underwent a series of remarkable outbursts in early 2002. 
During the period covered by our observations, 2002 December-2003 December, the 
near-infrared spectrum continued to show many of the characteristics of a very cool 
supergiant. However, throughout this period the spectrum also revealed strong and 
variable SiO first overtone emission, and Pa/3 emission. The 2003 December spectrum 
contained a series of Ti I lines with inverse P Cygni profiles. This is clear evidence that 
some material is falling inward towards the star. 
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1 INTRODUCTION 
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Iwisniewski et al.l 





Baneriee fc Ashokl |20j 

Rushton et al.l 120051) . The object was first detected in 
outbu rst on 2002 January 6 (V ma x = 10) ferown et alJ 
2002), and subsequent outbursts developed in 2002 February 
(Vmax = 6.7) and 2002 March (V m » = 7). Optical spec- 
troscopy showed a cool, reddened object, whose continuum 
shape and absorption spectrum were broadly consistent 
with those of a K supergiant. However, many of the spectral 
lines displayed P Cygni profiles, indicating outflow speeds 
of ~200km s _1 , and exhibited strong variability. The post- 
eruption phase was characterised by a dramatic increase in 
V — I, accompanied by a rapid trend towards later spectral 
types. As of 2002 May the visual magnitude had faded 
by ~ 7 mag since the 2002 Marc h peak, and the eff ective 
temperature was 2400 K iBaneriee fc Ashokl f2002). By 
2002 September further bolometric fadi ng had revealed the 
prese nce of a possible B3 V companion ( Dcsidcr a fc Munaril 
2002). At this time the near-infrared spectrum of V838 Mon 
displayed very deep H 2 0, CO, TiO, VO and AlO bands, 
and superficially resembled L -type brown dwarf s, indicating 
a very cool giant/supergiant fevans et alJl20o"3) . 

V838 Mon may not be a unique object. The similarities 
with V4332 Sgr and M31 RV are broad-based iMunari et alJ 
l2002l:lKimeswenger et al.l2002l ; lB"oschi fc Munarl2004l) . and 
a further ana logue may be found in the "Peculiar Vari- 
able i n Crux" JPella Valle. Saviane fc Wenderothl20o3) . Re- 
centlv lBaneriee^^dTT2^03i ~Dresented a near-infrared spec- 



trum of V4332 Sgr, which also displayed strong AlO bands, 
but in emission. Accretion models invol ving stellar merg- 
ers JSoker fc Tvlendall2003i) and planets iRetter fc Maroml 
2003) have been proposed to collectively explain this group 
of unusual objects. 

Here we present some of our latest observations of the 
near-infrared monitoring program of V838 Mon, emphasiz- 
ing the discovery of SiO overtone emission and inverse P- 
Cygni profiles on some atomic lines. 



2 OBSERVATIONS 

Low and high resolution infrared spectroscopy of V838 
Mon were obtained on numerous occasions between late 
2002 and end 2003, u sing the facility instrument CGS4 
jMountain et alJ ll99ol) on the United Kingdom Infrared 
Telescope (UK1RT). An observing log for the specific ob- 
servations presented here from that time interval is pro- 
vided in Table additional spectra obtained during that 
period will be reported later. The observations employed 
the CGS4 40 1/mm and echelle gratings at resolving pow- 
ers of ~1400 and 20,000 respectively. Flux calibration and 
removal of telluric absorption features were achieved by di- 
viding the target spectra by spectra of nearby calibration 
stars. To avoid spurious emission lines in the ratioed spec- 
tra due to hydrogen absorption lines in the calibration stars, 
the Pa/3, Pfy and Brcv lines were removed from the calibra- 
tion stars by interpolation prior to ratioing. The low reso- 
lution ratioed spectra were then multiplied by a blackbody 
corresponding to the effective temperature and broadband 
photometric magnitudes of the standard star to give the fi- 
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Table 1. Observing log. 



Observation date 


Wavelength coverage 


Resolving Power 


Std. Star 


(UT) 


(H 






2002 December 17 


2.9-4.1 


1400 


BS2530 (F2 V) 


2003 April 4 


2.9-4.1 


1400 


BS2530 (F2 V) 


2003 October 6 


2.9-4.1 


1400 


BS2530 (F2 V) 


2003 December 6 


1.278-1.286 


20,000 


BS1879 (08 III) 


2003 December 6 


3.989-4.021 


20,000 


BS1879 (08 III) 


2003 December 17 


4.016-4.046 


20,000 


BS1879 (08 III) 



nal flux-calibrated spectra. The high resolution spectra were 
scaled to make the continuum equal to unity after division. 
At low resolution, wavelength calibration was carried out by 
obtaining argon lamp spectra immediately before the target 
and is accurate to better than 0.001 (im. Argon lines also 
were used to calibrate the echelle spectrum at 1.28 /im, but 
telluric N2O absorption lines in the spectrum of the calibra- 
tion star were used to calibrate the 4 /im echelle spectra, in 
both cases to an accuracy of ±3 km s _1 . 



3 THE OH AND SiO BANDS 

The spectral evolution of V838 Mon (2002 December-2003 
October) in the range 2.9-4.1 /im is shown in Figure The 
bulk of the detail in the 2.9-3.5 /an region is due to a combi- 
nation of the Vi and W3 stretching bands of H2O. We have ob- 
served strong water bands in V838 Mon p r eviou sly in JHK 
spectroscopy, as discussed in IE vans et alJ i2003l) . 

The OH fundamental (Aw = 1) P branch lines strad- 
dle this wavelength range in late-type, O-rich stars. The 
line positions of the v = 1 — » band, obtained from 
iMaillard. Chauville fc Mantzl il976T) are indicated in the 
Figure. The line strengths are usually slowly varying across 
the band, as the excitation energy varies slowly as a func- 
tion of J. The P(19) - P(23) lines are the only OH lines 
present in the spectra and their intensities clearly fade over 
the period of our observations. The lack of higher excitation 
OH emission lines can be ascribed to temperature, but it is 
difficult to account for features at lower J unless OH ab- 
sorption, increasing towards lower J, is also present. Higher 
resolution observations might reveal such absorption. OH 
transitions involving higher vibrational states do not appear 
to contribute significantly in this wavelength range. 

The SiO first overtone (Aw = 2) is conspi cuous in the 
A > 4 /im region in emission, as was reported bv lLvnch et alJ 
(2004). In all of our spectra the v = 2 — > band is clearly 
visible, but there is little or no evidence for the w = 3 — » 1 
bandhead. Therefore only the lowest vibrational states w < 2 
are significantly populated. 

To produce a significant population in the v = 2 level 
by collisional excitation from v = 1 the density must be 
>Aio/Ci2, where A10 is the Einstein coefficient for spon- 
taneous decay from w = 1 to w = 0. For SiO-H2 collisions 
we can estimate the rate constant C12 for collisional excita- 
tion from the v = 1 to the w = 2 level using the vib r ationa l 
relaxation time form ula given in iMillikan fc White! Il963t) . 
Assuming, as is usual JScoville et alJl 980). the de-excitation 
cross sections scale like the radiative matrix elements, and 
using detailed balance, we obtain C12 — 5.6 x 10 -13 cm -3 . 



The required density is then riH 2 ~ 9.2 x 10 12 cm -3 at the 
excitation temperature of the SiO (see below). This will be 
an overestimate if there is a significant fractional amount 
of atomic hydrogen (riH <^nH 2 /100), since laboratory exper- 
iments with CO show that the rate coefficents for collisions 
with H are ~ 100 t imes larger than those involving H2 (see 
IScoville et aljll980l and references therein) . 

Figure [5] shows a later spectrum at high resolution in 
the spectral range 3.98-4.05 /im. This covers portions of 
the v = 2 — > band of the main isotopomers 28 Si ls O and 
29 Si 16 0. The fine structure is now resolved and emission in 
the individual rotational-vibrational lines is clearly appar- 
ent; all of the SiO lines in this range are R branch transitions. 
The 28 SiO line positions used in Figure [5] were obtained from 
lLovas. Maki fc Olsonl il98lfl : these include the positions of 
the R banch lines extending from R(19) at 4.0452 /im to 
the bandhead at i?(68) (4.0043 /im). The locations of lines 
from higher J values are not shown; the wavelengths of these 
overlap those of the ascending R branch but the lines are 
undoubtably much weaker. The radial velocity of the low J 
lines is +79 ± 4km s" 1 (heliocentric). 

In Figure |2] we present a preliminary attempt to 
model the SiO emission, using a simple, optically thin and 
isothermal slab with a Boltzmann population. The tran- 
sition probabilities of the 28 SiO lines used in the model 
were calculated using the di pole matrix elements given in 
iTipping fc Chackerianl dl98ll) . The less-abundant 29 SiO iso- 
topic species is not included in the model. However, there 
is no clear evidence for its detection in this complex spec- 
trum. The 29 SiO bandhead occurs at 4.029 /im, and the 
terrestrial abundance ratio 28 Si/ 29 Si is ~ 20. The lack of 
isolated lines in the spectrum makes it difficult to yield a 
definitive isotopic ratio, and it is only reasonable to con- 
clude 28 Si/ 29 Si > 5 from the present data. 

We find that a model with SiO excitation temperatures 
of T ~ 1200 K is the best representation of the observed 
spectrum (a more detailed analysis will be presented else- 
where). This is consistent with the low temperatures seen 
elsewhere in the infrared spectrum jEvans et aLl feoOS). and 
suggests the SiO must be located close to the stellar photo- 
sphere, rather than further out in the extended atmosphere. 
However , a lower SiO excit ation temperature of 790 K was 
found bv lLvnch et alJ i2004|) from their earlier spectra (2003 
January and February). Our low resolution data clearly show 
the SiO emission varying over the twelve months to 2003 
December (see Figure |5J, with the strongest emission in the 
2003 April spectrum, and it is just as pronounced in the 
Lynch et al. spectra, albeit at a lower resolution (see their 
Figure 8). A comparison of their results with our 2002 De- 
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Figure 1. LL' band spectra of V838 Mon on the dates indicated. The vertical lines denote the positions of the OH fundamental 
transitions shown with their corresponding J values. The wavelengths of the SiO v = 2 — * 0, 3 — * 1 and 4-^2 bandheads are also shown. 
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Figure 2. (a) 2003 December high resolution spectrum of V838 Mon showing the SiO first overtone emission. The figure also shows five 
unidentified emission lines . (b) Best fitting model spectrum computed with an excitation temperature of T = 1200 K, superposed on 
the observed spectrum. The positions of the 28 SiO V = 2 — > Ft branch lines, redshifted so they occur at +79 km s — 1 (heliocentric), are 
indicated in each panel (see text for details). 
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Figure 3. A high resolution spectrum of V838 Mon in the region 
of Pa/3 on 2003 December 6. 

cember 17 spectrum, suggests that the flux in the v — 2 — > 
band had increased by ~ 2.5 times during the following three 
weeks. 

We note that our model spectrum does not match the 
observed spectrum at a number of wavelengths in the SiO 
band. This is almost certainly due to the presence of other, 
as yet unidentified, lines. Five strong unidentified emission 
lines can be seen in the region shortward of the SiO band 
head; their observed wavelengths are 3.9907, 3.9926, 4.0003, 
4.0026 and 4.0048 /im. It is probable that other such lines 
overlap the SiO band. 



4 EVIDENCE OF INFALL 

Figure [3] shows the high resolution spectrum obtained in 
the wavelength interval 1.278-1.286 ycm. Clearly visible are 
a sequence of three spectral lines bearing well-defined in- 
verse P Cygni profiles with their absorptions centered at 
1.2829, 1.2838 and 1.2854 fjxa. The only spectral features 
which can possibly account for these lines are the permitted 
transitions of Til at 1.282512, 1.283495 and 1.285054 /im. 
These resolved Ti I lines are all remarkably narrow, with an 
average deconvolved FWHM of 27 ± 4 km s" 1 . The cores 
of the absorptions occur at a heliocentric radial velocity of 
+88 ± 7 km s _1 ; the peaks of the emission components are 
located at +65 ± 6km s _1 . All of these Til lin es originate 
from the b3-z3fo excited states Ivan Hoodl2004h . The exci- 
tation energies of the b3 lower levels lie between 1.43 and 
1.46 eV above the ground state, and thus could be collision- 
ally excited in a low temperature gas. There is a further Ti I 
line (1.28115 (im) present in this spectral region at an ob- 
served wavelength of 1.28185 /im, but the lower level of this 
a3D-y3fo transition is further above the ground state (2.16 
eV). This line is in absorption, but may also possess a weak 
emission component. It is difficult to identify the continuum 
level in this spectral region. 



5 HYDROGEN EMISSION 

A notable feature of several low resolution J band spectra of 
V838 Mon that we have obtained has been the prese nce of 
an emission feature near 1.28 (jm jEvans et alJl2005t . This 



feature is present in the high resolution spectrum shown 
in Figure [3] The most logical identification for it is Pa/3 
(1.282159 fj,m). This identification has not been confirmed 
by detection of other lines. Paa (1.875613 /im) occurs in 
a region of heavy H2O blanketing in the star and strong 
attenuation by telluric H2O. Pa7 was also not seen in the low 
resolution data and was not expected given the weakness of 
the Pa/3 emission. Bra, another expectedly weaker hydrogen 
line, is not apparent in the low resolution data presented here 
in Figure 

The Pa/3 line peaks at a heliocentric radial velocity of 
+56 ± 9 km s _1 . The red wing of the emission is almost 
abse nt. This is a signature of an a ccelerating/ decelerating 
flow dBohm-Vitense fc Quercilll987l) . but the superposition 
of other spectral features could also account for the asym- 
metric profile. Pa/3 absorption, redshifted by +96 + 6 km s _1 
(heliocentric), is a possible explanation. 



6 DISCUSSION 

We have presented infrared observations of V838 Mon which 
show SiO overtone emission, Pa/3 emission and Ti I lines with 
inverse P Cygni profiles. 

There have not been many reported occurrences of 
SiO overtone emission in astronomical objects. Of the 
few known inst ances, the phenom enon has been observed 
in S N 1987A jMeikle et alJll993ft . the Mira variable, o 
Ceti l|Ya mam^ae^iT IT999l) . the RV Tauri star, R Set 
iMatsuur^^^ani2002TT ! and is suspected in the latest M gi- 
ant stars. 

The observed SiO bandheads and band structure of the 
M giant stars SW Vir and RX Boo have been shown to be 
too weak as compared with the predictions by model atmo- 
spher es with appropriate effective temper atures <|Tsuj i et all 
Il994h . Si milar results were noted by iRinsland fc Wina 
(1982) and lTsuii et alJ dl997l) . who suggest that SiO emis- 
sion occurs in the outer atmospheres of the coolest O-rich 
giants, and fills in the strong photospheric SiO absorption 
lines. 

The presence of inverse P Cygni profiles in the post- 
eruption spectra of V838 Mon is clear evidence that some of 
the expelled material is now falling toward the star. The 
presence of SiO and Pa/3 emission can be understood if 
this infall is being compressed and shock-heated by material 
close to the star, where there are the high densities neces- 
sary to collisionally excite the SiO (see §3). This is similar 
to the phenomenon seen in pulsating variables where rising 
and falling layers lead to repetitive shock and line emis- 
sion very close to the photosphere l|Matsuura et all 120021 : 
iFerlet & GillelJll984t Icillet et al.lll9~89h . 

The SiO emission is redshifted by the same velocity as 
the Til absorption and, if present, the Pa/3 absorption. We 
therefore conclude that these atomic features must (like the 
SiO) arise very close to the star. The velocity shifts displayed 
by the Til and Pa/3 emissions are similar. The simplest ex- 
planation is that Ti I and Pa/3 trace the same material, with 
a spherically symmetric distribution, but the difference in 
excitation suggests that Pa/3 arises in a deeper layer. Our 
present observations do not allow us to determine whether 
additional velocity components are present in the SiO emis- 
sion. 
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Knowledge of the stellar velocity is required in order 
to convert our measured velocities to the rest frame of the 
star. This parameter is uncertain due to the presence of 
spectral lines with P Cygni profiles in the outburst spec- 
tra, and the non-det ection of mol ecular rotational em ission 
jRiishton et al.l2003h . Nonetheless iKipper et ail (|2004h sug- 
gested that the systemic velocity could be +59 ± 6 km s _1 
on the basis of two such velocity components in the P Cygni 
emissions in the outburst spectra, //this is the stellar veloc- 
ity then the infall velocity is ~ 15 km s , but again we em- 
phasize the uncertainty in the systemic velocity. The many 
different outflow velocities and components displayed by the 
spectral features in the outburst spectra demonstrated the 
compexities of the gas motions in V838 Mon at that time. 

Titanium is the only atomic species with transitions 
involving low-lying states in the high resolution spectral 
range presented here. Follow-up observations are necessary 
to search for infall signatures in other regions of the spec- 
trum, and to monitor the effect of this infall on the future 
behaviour of V838 Mon. 
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